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Abstract 

The Hamiltonian formulation of a plasma four-field fluid model that describes col- 
lisionless reconnection is presented. The formulation is noncanonical with a corre- 
sponding Lie-Poisson bracket. The bracket is used to obtain new independent fami- 
lies of invariants, so-called Casimir invariants, three of which are directly related to 
Lagrangian invariants of the system. The Casimirs are used to obtain a variational 
principle for equilibrium equations that generalize the Grad-Shafranov equation to in- 
clude flow. Dipole and homogeneous equilibria are constructed. The linear dynamics 
of the latter is treated in detail in a Hamiltonian context: canonically conjugate vari- 
ables are obtained; the dispersion relation is analyzed and exact thresholds for spectral 
stability are obtained; the canonical transformation to normal form is described; an 
unambiguous definition of negative energy modes is given; and thresholds sufficient for 
energy-Casimir stability are obtained. The Hamiltonian formulation also is used to 
obtain an expression for the collisionless conductivity and it is further used to describe 
the linear growth and nonlinear saturation of the collisionless tearing mode. 
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1 Introduction 



Due to the lower dimensionality of configuration space as compared to phase space, fluid 
models of the plasma have an intrinsic computational advantage over kinetic models. For 
this reason it is of great interest to develop and apply such models even in cases where the 
collisionality is too small to provide a firm justification for their use. In particular, fluid 
models have made important contributions to the understanding of magnetic reconnection, 
a phenomenon that plays a key role in events such as solar flares, magnetospheric substorms, 
and sawtooth oscillations in tokamaks. [U, Ej [3[ llj They have also made key contributions 
to the understanding of plasma turbulence in the core O [6], 13 E] and edge[9|, [TOl [TTl |12] 
of magnetic confinement experiments in coUisionless as well as collisional regimes. More 
generally, they offer the promise of being able to perform simulations of multiscale phenomena 
that are beyond the reach of kinetic models even after accounting for foreseeable advances 
in computation speed. [13] 

An important consideration when constructing new plasma fluid models is the existence 
of a Hamiltonian structure (see [lll[T5l[T6] for reviews). Because the fundamental laws gov- 
erning charged particle dynamics are Hamiltonian, dissipative terms, which ultimately arise 
from simplifications, must be accompanied by phenomenological constants such as resistivity 
and viscosity. When such phenomenological quantities are neglected, it is desirable that the 
resulting model be Hamiltonian, as is the case for the most important kinetic and fluid models 
of plasma physics. The preservation of the Hamiltonian structure provides some confidence 
that the truncations that are used to derive the fluid model have not introduced unphysical 
sources of dissipation. The presence of the Hamiltonian structure has the additional ben- 
efit of providing important tools for calculations. For example, the magnetohydrodynamic 
(MHD) energy principle is a consequence of the Hamiltonian nature of MHD. 

Since the discovery of the noncanonical Hamiltonian structure of MHD [17] , many plasma 
fluid models have been shown to possess a Hamiltonian description in terms of noncanonical 
Poisson brackets, e.g. [ll[6l[7l[8l[T8l[T9l[20l[2ll[22l[23]. In some cases, the requirement 
that the dynamics has Hamiltonian form has been used to guide the construction [21 [20] and 
has led to the identification of new and physically important terms. [20] In another case, the 
absence of a Hamiltonian structure for a given model was shown to lead to the violation 
of the solubility conditions for the equilibrium equations. |24j The Hamiltonian structure of 
fluid models has also been shown to be important for the consistent calculation of zonal flow 
dynamics. [6], O [8] 

Several fluid models have been proposed to study electromagnetic plasma dynamics 
(see [251 ES] for reviews). Some of these models have been instrumental in advanc- 
ing our understanding of magnetic reconnection. [201 [251 EH H] In particular, the model 
of Ref.[20] led to the discovery of fast (compared to Sweet-Parker) magnetic reconnection 
by Aydemir. [28l [29] This model included the effects of finite ion temperature, but neglected 
electron inertia. An alternative model that neglected ion temperature but did include elec- 
tron inertia as well as curvature effects was proposed around the same time by Hazeltine 
and Meiss (HM). The HM model was originally used to provide a unified description of the 
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formation of current channels in semi-collisional and coUisionless regimes. Its Hamiltonian 
nature, however, has not been investigated until now. 

Interest in the effects of electron inertia was recently revived by a controversy over its 
influence on the rate of coUisionless magnetic reconnection (CMR). In CMR, the "frozen- in" 
condition of MHD is broken by the inclusion of electron inertia instead of resistivity. [301 [31] 
This led Schep and collaborators to study two models that may be viewed as limiting cases 
of the HM model (aside from the fact that they avoid the Boussinesq approximation, re- 
taining instead the density n in the form logn/no).[ll [211 l2Z] These authors constructed a 
Hamiltonian formulation for these models, which were subsequently used to demonstrate the 
role of phase mixing of the Lagrangian invariants during fast reconnection [22] and to inves- 
tigate the role of instabilities of the nonlinearly developed current sheet. [32] More recently, 
Fitzpatrick and Porcelli (FP) have considered another limiting form of the HM model that, 
compared to the model derived by Schep et al, is valid for a wider range of values of (3, the 
latter indicating the ratio between the plasma pressure and the magnetic pressure based on 
the toroidal guide field. The FP model also extends the models of Schep et al. by including 
the effects of parallel ion compressibility. It has subsequently been used to study two-fiuid 
effects on the Taylor problem|^33] and on the linear growth of tearing modes. [34j Rogers et 
al. have recently shown that the predictions of the FP model for the linear growth rate of 
the tearing mode are in good agreement with those obtained with the gyrokinetic code GS2 
when the ion temperature is not too large. [35] 

In the present paper, we investigate the Hamiltonian structure of the FP version of 
the model of HM.[25l [36] The paper is organized as follows. In Sec. [3] we present the 
noncanonical Poisson bracket and show that this bracket produces the equations of motion 
with the appropriate Hamiltonian. In Sec. [H we use the noncanonical Poisson bracket to 
obtain four infinite families of new Casimir invariants, three of which suggest that specific 
combinations of the field variables are Lagrangian invariants. In terms of these variables, 
the equations of motion and the Hamiltonian structure achieve a much simplified form. A 
preliminary version of these results was announced in [37] . 

The remaining sections describe a variety of applications that rely on the Hamiltonian 
structure. In Sec. [5] we describe a variational principle for equilibria of the system and 
show that they are governed by a generalized Grad-Shafranov system of a pair of coupled 
elliptic equations. We treat two examples: dipole equilibria with Bessel function solutions 
and homogeneous equilibria that support wave motion. Sec. [6] explores the latter example 
in further detail by showing how to construct conventional canonical variables for the linear 
dynamics. We show that the system possesses Alfven-like and drift-shear modes, obtain 
exact stability thresholds, and give a definition of negative energy modes. We note that the 
Hamiltonian form is indispensable for an unambiguous definition of negative energy modes. 
We also obtain energy-stability conditions, sufficient conditions for stability akin to the 6W 
criterion of MHD. Section [7] contains a derivation of the coUisionless conductivity that relies 
on the Jacobi identity of the Hamiltonian formulation, which we then use to obtain the 
tearing-layer parameter A' and the growth rate for the coUisionless tearing mode. In Sec. [8] 
we use the conservation of a Casimir invariant to obtain the nonlinearly saturated current 
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profile and compare it to that obtained by Rutherford. [55] In Sec. IH] we summarize and 
conclude. 



2 Model equations 



The model of [36] is given by the following equations: 



(1) 




(2) 



(3) 



(4) 



Equation ([1]) is a reduced Ohm's law where the presence of finite electron inertia, which makes 
it possible for MR to take place, is indicated by the terms proportional to the electron skin 
depth de- Equations ([2]), and (jl]) are obtained from the electron vorticity equation, the 
vorticity equation, and the parallel momentum equation, respectively. 

Considering a Cartesian coordinate system {x, y, z) and taking z as an ignorable coordi- 
nate, the fields Z, ip and v are related to the magnetic field B and to the velocity field v by 
the relations B = Wif) x z + (5*^°^ + cpZ)'!, and v = —SJtp x z + f z, respectively. Here B^^^ is a 
constant guide field, whereas cp = (3/{l + (3) and dp = diCp with di indicating the ion skin 
depth. For small /?, dp ^ ps, the sonic Larmor radius. The ions are assumed to be cold, but 
electron pressure perturbations are taken into account and are given 

with P^'^^ a constant background pressure and pi coupled to the magnetic field via the rela- 
tion pi ~ —cpZ. Notice, here the parameter /? is defined as /? = (5/3)P^°-'/-B^'^^ , and above 
all the quantities are expressed in a dimensionless form according to the following normal- 
ization: V = aV, t = VAt/a, B = B/Bp, where a is a typical scale length of the problem. Bp 
is a reference value for the poloidal magnetic field, and va_ is the Alfven speed based on Bp 
and on the constant density. Finally, [f,g] '■= V/ x Vg ■ z, for generic fields / and g. 

3 Hamiltonian formulation 

A desirable property for fluid models of the plasma is that the non-dissipative part of their 
equation of motion should admit a noncanonical Hamiltonian formulation [HI [151 [IS]- 
short this means that it is possible to reformulate the ideal part of an n-field model as 



where C,i are suitable field variables, H is the Hamiltonian functional and {, } is the Poisson 
bracket consisting of an antisymmetric bilinear form satisfying the Jacobi identity. 



dt 



I = 



(5) 
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The first task in the derivation of a noncanonical Hamiltonian formulation is to iden- 
tify a conserved functional, usually the energy, that can serve as the Hamiltonian of the 
model. If one considers, for instance, a square domain T) in the x — y plane with doubly 
periodic boundary conditions, the four-field model ([l])-(!4]) admits the following constant of 
the motion: 

H = ]- I d?x {dlJ^ + I V^p + I V<^p + v^ + Z^) (6) 

with J = — V^V indicating the parallel current density. The quantity H represents the total 
energy of the system. The first term refers to the kinetic energy due to the relative motion of 
the electrons with respect to ions along the z direction. The third and fourth terms account 
for the kinetic energy, whereas the second and last terms account for the magnetic energy. 

Adopting ipe = '^ — dl'V'^'ip, U = V^<^, Z, and v as field variables, i.e. ^ = {ipe, U, Z, v), and 
([6]) as Hamiltonian, it is possible to show that the model can indeed be cast in a noncanonical 
Hamiltonian form with the following Lie-Poisson bracket: 

{F, G} = j d'x {U[F^, G^]u + MF^, G^]^^ + Z[F^, G^]z + v[F^, G^],) , (7) 

where 

[F(^,G^]u = [Fu,Gu] 

[F^,G^]z = [Fz,Gu] + [Fu,Gz] — di3de^[F^^,G^^] (8) 

+Cfsd,^{[F,, G^J + G,]) - a[Fz, Gz] - cp^[F,, G,] 
[Fg, G^]^^ = [F^^,Gu] + [Fu, G^J - dp{[Fz, G^J + G^]) + c^([F„ Gz] + [Fz, G,]) 
[F^, G^], = [F„ Gu] + [Fu, G,] + Cfsde\[Fz, G^J + [F^^,Gz] - cpj{[F,, Gz] + [Fz, G,] , 

with a = dp + cpde^/di, 7 = d^ /di, and subscripts indicate functional differentiation. It is 
straightforward to show that Eq. ([5]) with the Hamiltonian given in Eq. ([6]) and the above 
bracket reproduces equations of [36] • The Jacobi identity for the above Poisson 

bracket is readily established by applying the method described in Refs. [HI IS^ . 

Because the number of parameters is escalating, we record their definitions here for later 
referral: 

dp = cpdi , d = \J d1 + df ^ 'y = d^/di^ a = dp + cpdj^ /di = cpd'^ /di . (9) 

The basic parameters of the model are cp, de, and di, while dp, a, 7, and d are useful short- 
hands. 



4 Casimir invariants and bracket normal form 

Lie-Poisson brackets for noncanonical Hamiltonian systems are accompanied by the presence 
of Casimir invariants. A Casimir invariant is a functional that annihilates the Lie-Poisson 
bracket when paired with any other functional. That is, a Casimir G must satisfy 

{F,G} = 0, (10) 
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for every functional F. Thus, Casimir invariants constrain the nonhnear dynamics generated 
by the Poisson bracket for any choice of Hamiltonian. 

In order to identify the Casimirs of the four-field model we proceed in the following way. 
First, multiplying Eq. (jlj) times di and adding it to Eq. ([T]), we find 

dD 

-g^ + [v,D] = 0, (11) 

where D = ijje + diV is the ion canonical momentum. Equation ( ITTi) indicates that the field 
D is a Lagrangian invariant that is advected by the flow generated by the stream-function ip. 
The presence of this Lagrangian invariant also suggests that using D as one of the variables 
will simplify the Lie- Poisson bracket. Indeed, upon replacing ip^. with D as field variable, 
Eq. (ITU]) for the four-field model becomes 

{F, C} = j d\ {Fu[Cu, U] + Fd[Cu, D] + Fu[Cn, D] 

+cpF,[Cz, D] + cpFz[C,, D] + Fz[Cu, Z] + Fu[Cz, Z] 

-aFz[Cz, Z] - Cf,^F,[C,, Z] + F,[Cu, v] + Fu[C,, v] 

-aF4Cz,v] - aFz[C,,v]) = 0, (12) 

a simpler bracket. 

For Eq. ffTOj) to be satisfied for any F, it is necessary for the coefficients of each of the 
functional derivatives of F in f[T2]) to vanish separately. This leads to the following system 
of equations for C: 

[Cu,D] = 0, (13) 
[Cu,U] + [Cn,D] + [Cz,Z] + [C,,v] = 0, (14) 
-cp[C,,D]-[Cu,Z] + a{[Cz,Z] + [C,,v]) = 0, (15) 
C/,[Cz,D]-cp^[C,,Z] + [Cu,v]-a[Cz,v] = 0. (16) 

The problem of finding the complete set of Casimir invariants is thus equivalent to solving 
the set of four coupled equations (fT3|) - (fT6|) . 

Beginning with (fT3|l . a functional integration shows that C is of the form 

C{U, D, Z,v) = j d^x {UT{D) + g{D, Z, v)) , (17) 

where JF and g represent arbitrary functions of their arguments. The problem is now reduced 
to finding the functions T and g. Considering next Eq. 0141) . we see that this equation is 
automatically satisfied for any choice of C with an integrand that depends only upon the field 
variables and not their spatial derivatives, and therefore imposes no constraints. Substitution 
of ([nD into ([ISD yields 

{cfig^y + ag^D) [f , D] + {cjsgvz - T'{D) + agoz) [Z, £>] = , (18) 
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where ' indicates derivative with respect to argument and subscripts on g denote partial 
derivatives. In f|T8|) the coefficients multiplying the brackets ' [ , ] ' must vanish independently. 
These two relations lead to 

C(3g. + ago = ZT'{D) + K{D), (19) 

where K an arbitrary function of and integration of (|T9l) by the method of characteristics 
gives 

g{Z, V, D) = -J'iD) + JCiD) + g{Z, v - cpD/a) , (20) 
a 

where /C' = K/a. Similarly, insertion of f[T7l) into f[T^ yields 
{cpgzv - J^'{D) + agzo) [D, v] + cp {gzz + 19vd) [D, Z] + (c^75'„„ - agzz) Z] = , (21) 

which gives three equations, but only gzz + IdvD = provides new information. Defining 
V = V — c^D/a, we see that Q satisfies the wave equation 

Qzz-''-^Qvv = ^. (22) 
a 



and thus has the general solution 



^ V 2a V 7« V J 



^]-^^■ (23) 



2a J 



Therefore, in light of ( fTTl) . ( !20l) . and ( |23l) we have the following four independent families 
of Casimir invariants: 

Ci = [ £x (u + -] J^{D) , (24) 



Ca = y SxK.[D), (25) 

where we have scaled the arguments of C± for a reason that will become apparent soon. 

Knowledge of the functional dependence of the Casimirs suggests a simplification of the 
Lie-Poisson bracket will occur if it is written in terms of the new variables 

D = D, (27) 

C = U + -, (28) 

a 

T± = ±Tr4sX i^i^ - T ^^-^) ' (29) 
zc/sd-^df. ^ ' 
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which possess the inverse relations 



D = D, (30) 
U = ( + T+ + T^, (31) 
Z = -a{T+ + T_), (32) 

V = ^D-^(r+-r_). (33) 

Indeed, in the new variables, the Lie-Poisson bracket reads 

{F, G} = J d'x iC[Fo Gd + Di[Fn, G^] + [F^, On]) 

+T4Ft_,Gt_]+T4Ft^,GtJ) , (34) 

which is a bracket normal form that relies on the scaling used above. The bracket in terms 
of the new variables reveals its algebraic structure: it is identified as a sum of direct product 
and semi-direct product parts [191 [IHl ES] and, consequently, the Jacobi identity follows from 
general theory. Making use of the variables suggested by the form of the Casimirs, the model 
equations can be rewritten in the compact form 

dD 

^ -b,^], (35) 



dt 
dC 
dt 

where for convenience we have defined 



-[^X]+d-'[D,i;], (36) 
-b±,T±], (37) 



:= ± . (38) 

Note, the variable ( plays the role of a "generalized" vorticity, and our development reveals 
the existence of the three Lagrangian invariants D, T+ and T_ associated with the families 
of Casimirs G2, C+ and C_, respectively. The existence of such invariants implies that the 
values of D, T+ and T_ are constant on the contour lines of ip, ip+, and respectively. By 
choosing "top-hat" functions for the free functions in the Casimirs G2, C+ and C_, it follows 
that the area enclosed by the contour lines of the Lagrangian invariants remains constant. 
Notice also that T+ and T_ in the limit j3 —>■ and di ^ 00 become proportional to the 
Lagrangian invariants G± = ip — d^'V'^tp ± dePsU of the two-field model derived in [4]. The 
family Gi is of a different nature and one of the constraints imposed by it is that the total 
value of ( within an area enclosed by a contour line of D remains constant. 



5 Equilibria 

Equations governing the equilibrium of the FP model are most easily obtained by setting 
the time derivatives equal to zero in Eqs. (I35|) - (l37|) and solving for the the fields D, (, 
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and T±. Alternatively, it is possible to derive equilibrium equations from a variational 
principle, the existence of which is assured by the Hamiltonian nature of the equations. [15] 
The construction of the variational principle is more laborious than the direct approach but 
the extra work is richly rewarded by the well-known benefits of variational principles. In 
particular, the variational principle provides a basis for studying the stability as well as the 
equilibria of the system. 

For a system with a collection of Casimirs, extrema of the free energy functional F = 
H + are equihbria of the system. Such extrema can be derived by setting the first 
variation 5F equal to zero. If ^j, i = 1, 2, . . . , denotes the field variables of the system, this 
amounts to solving the equations H^^ + ^ C^. = 0, where the subscript indicates functional 
derivative with respect to ^i. One advantage of this variational approach is that for equilibria 
obtained as extrema of F, the second variation of F provides criteria sufficient for stability. 



5.1 General equilibria 

For the model here, the free energy functional reads 



d X 



+ 



2 
d^ 



y2 ^2 

y + — + /C(V'e + div) 



di 



di 



2cRdd. 



■V — 



2cad^ 



(39) 



where the linear operator C is defined by C^^ip = ip — d'^V'^ip = ipe- Notice that the 
arguments of the functions present in the Casimirs are of course much less compact when 
written in terms of the variables ^ = {ipe,U, Z,v). On the other hand, in terms of the 
variables T := (D,(^,T+,T_) the free energy functional reads 



F[T] 



d'^x 

V 

1 fd, 



c^d^ 

-^(t2 + t!) + 



D + cpd{T+ -T_)]C(-^D + cpd{T+ - T. 



(40) 



Equation (HOl) shows that in terms of the variables that are "natural" for the Casimirs, 
the expression for the Hamiltonian becomes complicated. Unfortunately, there exists no 
preferred set of variables in terms of which both the Hamiltonian and the Casimirs take a 
simple form. In order to obtain equilibrium solutions by means of the variational principle, 
it is convenient to choose the T variables that are natural for the Casimirs, calculate the 
required functional derivatives of the Casimirs with respect to these variables, and then use 
the functional chain rule to obtain the functional derivatives of H in terms of the variables 
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^. More specifically, by setting 5F = 0, the resulting equilibrium equations are given by 

= H^ + J2C,^ = H^ + HD)=0, (41) 

j 

Fd = HD + Y,C,o = HD + CJ"{D) + IC'iD) = 0, (42) 

j 

Ft^ = HT^ + J2c,^^=HT^+g'^iT^) = 0. (43) 
j 

where the index j ranges over the set {1,2,+,—}. The functional chain rule then can be 
used to evaluate the functional derivatives of the Hamiltonian, 

= Hu, 

Old 

Ht± = ±cpded H^^ + Hu - aHzT Hy. 

The functional derivatives H^^, Hu, Hz and H^ can themselves easily be obtained from the 
Hamiltonian written in the form ([6]). The equilibrium equations (l4T!) - (H3l) are then given by 

-y^ + HD) = 0, (44) 
-§^'i^ + ^ + CJ"{D) + }C'{D) = 0, (45) 

^cpdJV^-V-aZT^ + G'^T^) = 0. (46) 

These equations are expressed in a mixture of the ^ and T variables. In order to simplify 
them we eliminate V'^ip from the last two equations by using 



dy^i) = ip- D + diV. 
Using this in Eq. (H5|) - (H6|) . we find 

-d-^^ + C^\D)+t\D) = 0, (47) 
-i^±±^{diD-d\Tdd,Z) + g'^{T^) = 0. (48) 

where fC{D) = IC{D)+D^/2d^. Equation ( HHl) may be simplified further by expressing Z and 
V in terms of the T variables using (l32l) - (!33l) . This leads to the following complete system 
of equilibrium equations: 

-ip + J^{D) = 0, (49) 
-y^± + Q±{T^) = 0, (50) 
-d-^ij + Cr{D) + }C\D) = 0, (51) 
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where G±iT±) = G±{T±) +a'^T^. One easily verifies that Eqs. P9|) - fl5Tl) describe equihbrium 
states by substituting them into Eqs. fl5^ - fl57l) . 

Continuing our present approach of ehminating the ^ variables would now lead us to 
express the ip, ip, and (p± in terms of integral operators acting on the T variables. Clearly 
this is undesirable. Instead, we note that the above four equations express a dependency 
between the six quantities V^ip, V^V9, ip, ip, Z and v. It is thus possible in principle to use 
these equations to express four these quantities in terms of the remaining two. If we choose 
(p and ip as the independent fields we will obtain a closed system of equilibrium equations of 
the form 

= Siij,^), (52) 
VV = P{^,^), (53) 

Equation (l52l) is a generalized version of the Grad-Shafranov equation, whereas (l53i) is an 
analogous equation that determines the equilibrium polarization. 

In order to calculate the form of the functions S and P we invert Eq. (H9|) - (l50|) : 

D = a{ip), (54) 
T± = t±(<^±), (55) 
C = d-2^a'(<^) + %), (56) 

where a and t± are the inverses of and Q± respectively and b{ip) = —ft\a{ip))a\(p). Solving 
these equations for V'^ip and V^V then yields fl52l) - fl53l) with 

PW'.V) := %) + *+(»>+) (58) 

We complete the system by expressing the two remaining unknowns Z and v in terms of p} 
and Tp. From Eqs. fl32|) and fl33l) . using flM|) . fl55l) there follows 

v{il),p) = dia{p)/(f' - cpded{t+{p+) -t_{pJ))/di] (59) 
Z{ij,p) = -a{U{p+)+t_{p_)). (60) 

We have thus shown that solving the equilibrium equations amounts to solving the cou- 
pled system of fl52l) and (!53|) for the unknowns ip and p). This requires making choices for 
the free functions t±, a, and h. If one is only interested in solving the equilibrium problem, 
one may determine these free functions directly from physical considerations and the rela- 
tionship between these functions and those appearing in the variational principle may be 
ignored. These relationships become important, however, if one wishes to use the variational 
principle either to solve the equilibrium or the stability problem (the variational principle 
has well-known advantages both for numerical and analytic applications). In this case the 
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functions Q±-,T^ and AT, appearing in the variational form may be determined in terms of 
a, 6, and t± as described above. Variational treatments of two-fluid equilibria have been 
given by [IHl HI], |12] , and applications of variational principles to stability are discussed by 
[40l|4ll|43l|44|. 

We conclude this section by noting two difficulties with the equihbrium equations (l52i) - 
( l53l) . The first is that these equations may become hyperbolic in the presence of strong fiows. 
Recent analyses of this problem can be found in [121 135]. The second difficulty is that the 
right-hand sides are singular in the limits d^^L and di -C L, where L is a macroscopic scale 
length. [iO| HH US] That is, for macroscopic equilibria the derivatives in Eqs. (!52|) - (l53|) are 
multiplied by a small parameter, so that these equations form a stiff system. This has led 
to considerable grief, in particular for Field-Reversed Configuration (FRC) devices where 
a similar set of equations is encountered. Steinhauer has proposed a method for dealing 
with this problem that he named the "nearby fiuid" approximation. [l6l HT] In the following 
section we outline a similar approach to solving Eqs. fl52|) - fl53|) . 



5.2 Perturbative solution for macroscopic equilibria 

In order to deal with the singular nature of the equilibrium Eqs. (1521) - (!53l) . we expand the 
fields in powers of the small parameters de and di and solve term by term. A byproduct of 
this procedure is the clarification of the physical meaning of the profile functions. 
We begin by considering the limit rfe — > 0. For convenience we define 

i± i^i) + := deCf3dt± (v?±) 

and expand this and other profile functions in powers of de according to 

4(0) = (0) + d/l\<P) + diit\<p) + .... 

We then consider the equilibrium equations order by order in d^.. From fl53|) and fl58|) we 
obtain at lowest order = using this result in (152!) and (1571) we obtain to 

lowest order = £L°^(</)) = 0/2. To next order, ([52]) and ([57]) yield = 

substituting this result in (l53l) and ([58]) . we find the following equation for the vorticity 

VV = Kv) + a'{^)^/d,^ - h{tP)/dp , 

where 6 := 6 + y^/ci^ is a profile function for the vorticity and h{ilj) := —2t^'* (■?/'). 
From the terms of order unity in ( l52i) and ( 157]) we next find 

V^^ = -a{^)/d^ + h\i^)^/dp + /(V^), (61) 

where /(■?/') = t\_ {—ip) — t\ {ip). Equation ([6T]) is the Grad-Shafranov equation, where the 
term proportional to ip is the polarization current and I{ip) describes the inductive current. 
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The parallel velocity may be obtained from fl54p : to order (i^, D = ip—dlV'^tp+diV = a{(p), 
yields 

diV + ip = a{ip) . (62) 

Expanding fl60|) to order dg, gives 

Z + ^/dp = h{^) . (63) 

The sum Z+(p/d/3 represents the electron stream-function. The fact that the electron stream- 
function is constant on surfaces of constant flux, as expressed by Eq. fl63|) . is a statement of 
the frozen-in property. 

We may carry out the limit rfj ^ in a similar way. From the ion momentum conservation 
Eq. (162|) . we obtain ip = a{ip) showing that the electrostatic potential must be a flux function 
to lowest order. It is convenient to introduce := a~^{ilj) = ip to denote the inverse of 

a. We also define the Alfvenic Mach number M{ijj) := d^/dip. Note that Eq. (|63ll specifies 
that to lowest order, h^ip) = ^{%lj)/dp. In terms of these quantities, the vorticity equation 
shows that to lowest order. 

In order to eliminate (p from the Grad-Shafranov equation it is necessary to calculate the 
correction to the electrostatic potential. This is given by the vorticity equation, 

M'(VV')' + MV'V' = (M-2 - c^')^(2) _ 

Note that the potential exhibits a singularity for M = corresponding to the sound-wave 
resonance. Lastly, after eliminating the electrostatic potential from Eq. fl6Tl) we recover the 
MHD version of the Grad Shafranov equation, 

(1 - M^)V^^ - MM'{Vi)f = 

In the following sections we present some explicit solutions of the equilibrium equations for 
simple profile functions. 

5.3 Quadratic Casimirs— dipole equilibria 

The case of quadratic Casimir invariants is easily tractable. Choosing 

1C{D) = ^D\ T{D)=A^D, g^{T^) = ^Tl, (64) 
and following the steps of Sec. 15. II leads to 

«(.)^^. + M.J^^i^f^. (65) 
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Upon inserting fl^ into (jS^D and fl^ . we obtain 

V^tp = Sl^lJ + S2ip and V V = ^iV^ + ^2<^ , (66) 

where 5*1^2 and P12 are constants that depend on A±, A^, A^,, and the parameters of the 
system. Note, S'i_2 and Pi_2 are arbitrary except that 5*2 = —-Pi. Consequently, Eqs. ( l66l) 
have a variety of solutions that are closely related to the double-Beltrami flows investigated 
by Yoshida et al. [IHl HI] Specifically, Refs. [lOl HI] neglect electron inertial effects that are 
kept here, but they allow for more general geometry. In general, (!66|) can be diagonalized 
resulting in two decoupled equations of the form 

V'x. = -A.X«, ^ = 1,2, (67) 

where Ai,2 = — (•S'l + P2 ± a/ {Si — P2)^ — 4Pf )/2 and the Xi's are linear combinations of (p 
and ^z^. If a solution of this system is found, then one obtains v and Z as particular linear 
combinations of and ip as described in Sec. 15.11 

Rather than describe the general solution, we give an example representative of the 
kinds of solutions that are possible. We assume a circular domain of unit radius, adopt 
polar coordinates (r, 6*), and adjusting the parameters Si, Pi, and P2 so that a/Ai,2 are 
zeros (possibly distinct) of the first order Bessel function, i.e. Ji(a/Ai^) = 0. With these 
assumptions, we obtain the solution 

Xi{r,e) = A Jl{^/Xir) cose, (68) 

where the Aj's are constants and each of the x's has a dipolar structure like that depicted 
in Fig. m 

0.8- 
0.6- 
0.4- 
0.2- 
y 0: 

-0.2- 
-0.4: 
-0.6- 



-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 
X 

Figure 1: x contours for dipole solution of (1681) 
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5.4 Homogeneous equilibria 



The quadratic Casimirs of flMl) also yield homogeneous equilibria, i.e. equilibria for which 
the linear dynamics has constant coefficients. For this choice, the free energy functional of 
(140|) can be written as follows: 



V 



(69) 



where 
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Recall ^ = {ipe,U, Z,v) and T : 
where the matrix T is given by 



T 



and 



A 
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A+ 





^0 











(70) 



{D,C,T+,T^). Equations (127D-([29D amount to T = T^, 



/ 1 






1 








dj 
C0d^ 



di \ 




—dj dedj 
2c0d'^ 2ci3d^ J 



(71) 



The free energy functional can then be written as a quadratic form: 



V 



(72) 



whence the equilibrium equations are obtained upon setting the functional derivatives to 
zero, 

{H + r^AT)i = 0. (73) 

Here we have assumed the formal self-adjointness of the operators and V~^. Equation 
(173|l is a linear homogeneous system of four equations with the four unknowns ipe, U, Z, and 
V. The equilibria treated in Sec. 15.31 are solutions of this system, but of interest here are the 
homogeneous equilibria 



(74) 



where the is the Alfven speed, and az and describe density and velocity shear, respec- 
tively. These are clearly equilibrium solutions and can be related to the chosen Casimirs. 
Evidently, we are seeking solutions with V^ip = and f/ = 0, so (1731) reduces to 



(/° + r^iar)^ 







(75) 
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where ^3 := {iJ,Z,v), 



/O 0\ /Ad \ 

J° = 1 and is = A+ . (76) 

\0 1/ \ A_/ 

This is a hnear homogeneous system of three equations with unknowns ip, Z, and v. The 
existence of non-trivial solutions requires det(J2 + T^AT) = 0, which fixes a condition on 
Ad, y4+ and A- that assures ip, Z, and v are linearly dependent. Thus, they all can depend 
on X and be proportional, consistent with ([71]). These equilibria correspond to uniform 
poloidal magnetic fields, and to toroidal magnetic and velocity fields proportional to x. 



6 Normal forms for homogeneous equilibria 

Here we work out the linear canonical Hamiltonian form for the dynamics obtained by 
expansion about the equilibria of Sec. 15.41 In terms of the T variables the equilibrium is 

Do = ao X, Co = «c^5 ^±0 = a± X , (77) 

whence e.g. ao = 0.^ + dia^ and ( j33l) give = c^ddg^a^ — a-) + aoi^ — d"^ / d"^) . We obtain 
the Poisson bracket and Hamiltonian for the linear dynamics, obtain the dispersion relation, 
then we find a set of canonical variables, and discuss canonical transformations to normal 
forms. En route we define the meaning of negative energy modes. 



6.1 Linear Hamiltonian form 



Denoting the linear variables with a tilde, i.e. C = Co + C; stc, and expanding the bracket 
( l34l) gives the bracket for the linear dynamics (see [T5] ) 



d'x C( 



6F 5G' 
^C' ^C. 



5F SG 



+ 



5F 5G 



+ T. 



±0 



6F 6G 



where a sum over the ± terms is implied. Upon integration by parts this bracket becomes 



{F, G}i 



d X 



6F d 5G 



6C dy SC 



6F d 6G 5F d 5G 



ac_^^— + I — H --^ + a± — 



5D dy SC 6C dy 5D 



6F d 5G 



''5T+dyST. 



• (78) 



The above bracket together with the following quadratic Hamiltonian: 



Hl = ^ I d''x{ iJeJ + I V0p + v^ + Z^ + AdD^ + A±Tl 



where 



A 



D 



00 

\a£, 



— 1 1 and y4± = I ± 



d «v 2^^'^ 
de a± df 



(79) 



(80) 
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when written entirely in terms of the variables {(,D,T^,T_), yield the linearized equations 
of motion in Poisson bracket form. 

Unlike the nonlinear semi-direct product bracket of (1341) . the linear bracket of (ITS]) can 
be brought into direct product form by the transformation 

D = -a^i) + aDC , (81) 

which yields the convenient form 



{F,G}l = - j (fx 
where := — 



' 5F d 5G 6F d 6G 6F d 6G 



(82) 



6.2 Canonical coordinates 

Because the equilibrium equations do not depend explicitly on x and y, we Fourier expand 
C as 

oo 

ax,y,t)= Ck.,ky{t) e-^^'-^^'yy^ (83) 

and similarly for D, T+, and T_. For convenience we suppress the sum over and set ky = k. 
The variable k^ will only appear in the combination k"^ := k'^ + ky. It is not difficult to prove 
(see e.g. [48]) the following general functional derivative relationship: 

where F[(] = F[(k\. Note, in the above expression there is an extra factor of 27r that occurs 
in the denominator of the suppressed sum on k^, but this factor is compensated by a factor 
of 27r that accompanies a suppressed sum in the Hamiltonian. 

Upon inserting (18^ and similar relations for the other variables into the bracket of (178]) 
gives 

_ ^ik\ fdF dG dF dG\ ( 5F 5G 5F 5G\ 

dF dG dF dG \ ( dF dG dF dG 



where in order to facilitate the transformation to canonical variables the positive and negative 
values of k are assumed to be independent and the sum is written over only positive values 
oik. 
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A set of real valued canonical variables is given by 



with the inverse transformation 



_ 1 ka^ / (2) (2)\ ^ (2) 



1 ka+ ( (3) . (3)\ 1 / (3) , . (3) 



1 / (4) . (4)\ _ 1 ka_ / (4) (4) 



Above we have assumed that and a_ are positive; in light of its definition, is 

negative and there is an intrinsic parity difference between the D±k and the C,±k degrees of 
freedom. This is a fundamental property of linearized semi-direct product brackets. If a^, is 
positive, then is negative, and one must alter the C,±k transformations. In general, if any 
of the a's are negative, then a suitable canonizing transformation is obtained by inserting 
an absolute value inside the square root and replacing the corresponding by minus qk- 
The Hamiltonian corresponding to fl85|) is 

00 

0-1 . - 12 Oi\ 



Hl = 27r^4l^fcl'-— vl«sl/«c(^fcC-fc + 5_fea) + («2 + «i|«Bl/«c)ICfc|' 
k=i '^C '^C 

+ («2 + ^^I^sIk) + C-feT^,) . (88) 

In term of the canonical variables this Hamiltonian becomes 

fc=l ij'=l 



where the symmetric matrices M and K are given by 

Mil = Kii = -^^ai , M22 = K22 = ^ai + — , M33 = K33 = a+ag 



M44 = K44 = a_a4 , M12 = -A12 = 2~^i ' = --^13 = «6 



Mi4 = -f^M = - ^ ' " ae , M34 = ^^34 = ya^aTas , (90) 



with 



M23 = -ft'23 = 1/070702 ae , M24 = = Jaca_a2 H ae 



1 d\ 



"2 - T2- ' 



d? d\l + and? d\l + rf^^i ) ' 

1 

= ± . 4 24^^ _ 1 cpda^ 4d^ 

kl'^^^ dj 1 + C/2A;^ kl dea± 1 + Ci2fc2 
1 



05 = 7T + 



e 



= (1 + d{ki)d ■ ^''^ 

Note, the matrices M and have commuting and anti-commuting parts, which can be 
traced to the intrinsic parity difference mentioned above. 



6.3 Stability, signature, and normal forms 

One could proceed directly by linearizing the system of Eqs. about the equilibrium 

([71]) and obtain a system for the linear dynamics. Because we are using variables indexed 
by k and — as independent variables, we obtain the combined linear system 

ik = 'R-k- ik and = 'R-k ■ i-k , (92) 

where C,±k is a four dimensional vector and TZ±k is a 4x4 matrix, for each value of A; = 1, 2, ... . 
This means we have an eight-dimensional system for each value of k, and upon assuming 
^ ~ exp(ia;t) we obtain the dispersion relation from 

det I ^^"^^^ '^'^ . r^^rr, I = det(2Cij/4 — T^fc) ■ det(2Co'/4 — 7^_A:) 
\^ Oi iujU — rC-k J 

= det{-uj^h + nl) = 0. (93) 
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where O4 is a 4 x 4 matrix of zeros, I4 is the 4x4 identity matrix, and use has been made 
of TZ-k = —Tik which is easily verified for the case at hand. Alternatively, one can assume 
{Ok iVit) — (^fc^ 6xp(za;t),p[!^ exp(iti;t)), and obtain the dispersion relation from Hamilton's 
equations in the form 

lujqf = ^ kMij p^i^ and lupf = - ^ kK^^ q^^ , (94) 
i 3 

whence one obtains the dispersion relation as (\ei{uj'^Ii — k'^ M K) = 0, a relation equivalent to 
(l93ll with MK = TZ\. Because of this special form, the dispersion relation can be factored, 
and reduced to a 4 x 4 determinant that provides the frequencies for both positive and 
negative k. Thus the symmetry TZ^k = ~TZk allows for a simplification that is manifested 
by the special form of the Hamiltonian of fl89|) . We note that this special form occurs for all 
the basic fluid and plasma models, because they have real variable Hamiltonian form. 
For convenience, we introduce the dimensionless variables 

va = a^, K_L = k±de , r = — , M 



dp kvA 



_ de _ [rn^ _ aydpr _ azdpr 
di V i^i Va Va 



in terms of which we derive the following dispersion relation: 



k\)N^ - uN^ - {6^ + 6\l + 6s + r^ + 4) + J^r^N + r^6{s + 5) = , (96) 



From the form of (l96l) it is clear that is a function of alone; moreover, it can be 
rewritten as 



where 



Ar, = ^i±V^^!±M±3, N. = r, and Ns = Vl + S^ . (98) 

Analysis of the dispersion relation reveals that the four roots correspond to two that are 
Alfven-like, two that are a combination of a drift-like wave that arises from density shear u 
a Kelvin-Helmholtz-like wave that arises from the parallel velocity shear s. 

For example, if one sets S = and 1^ = 0, then (l97|l yields in dimensional variables 



uj = ±kvAy + k\d'^^)/{l + k'j_dl). The square root in this dispersion relation displays 

the slippage of flux, which comes from two sources, a numerator that is dependent on the 
electron temperature through dp and a denominator that is dependent on electron inertia 
through df.. Both of the terms in the square root break the MHD frozen-in condition and 
for de = dp = one recovers the Alfven wave dispersion relation. For de = and dp this 
dispersion relation reduces to that for a version of the 'kinetic' Alfven wave, while for large 
k one obtains the phase velocity a;//c ~ VAdp/de = \fT^fr(ie^ the electron thermal speed. 
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For = and large k, the lower hybrid frequency u ~ VA/de = VAOJpe/c = eB /cy/mjnl 
is obtained. Similar limits reveal the presence of drift-waves associated with u and Kelvin- 
Helmholtz modes associated with s. For example, if we set ^ = and suppose N is small for 
large wave-numbers, then (p7|l gives the drift-wave dispersion relation, uj = kv^/{l + k\d'j^), 
where = d^az- 

Equation ( |97j) is convenient for obtaining stability criteria, by examination of the zeros 
and divergences of its right hand side, and by noting that it asymptotes to unity for large 
|A^|. Because 5 and r are always positive, the sign of the divergence at = is governed 
by (5 + s, and this sign distinguishes two cases. The first case, 5 + s > 0, corresponds to 
positive or weakly negative parallel velocity shear. It is convenient to define the central band 
by C<5 := {A^||A^| < N^}, which is bordered by the divergences, and the set of frequencies 
at which zeros occur, M := {A^+, A^_, Nr, —Nr}. If any two elements of A/" are contained in 
Cs, then the system is stable. If three elements of J\f are contained in Cs, then the system is 
stable, and if all four of the elements of A/" are contained in Cs, then the system is stable. The 
second case, 6 + s < 0, corresponds to strong negative velocity shear. If u"^ + 46(6 + s) > 0, 
then the system is unstable for large enough k±, but always possesses two stable modes. In 
the case of very strong negative velocity shear u'^ + 6{6 + s) < 0, the set of zeros becomes 
M := {Nr, —Nr}, and the system is unstable for all n±, but again always possesses two 
stable modes. 

In Figs. [2] and [3] two solutions of (|97|) . with the real part of N versus k,±, are plotted. In 




0' + 2 4 6 8 10 



Figure 2: Solution of the dispersion relation of ( l96ll . The real part of is plotted vs. k± for 
r = 1.11, u = 10.5, S = 0.10, s = 18.5, which corresponds to positive parallel velocity shear 
with 5 + s > 0. The energy signatures of the modes are indicated. 
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Fig. [2] the parameters r = 1.11, 
the same values except s = — 13 



u = 10.5, S = 0.10 , and s = 18.5 are used, while in Fig. [3] 
.4 are used. For small k±, the uppermost curve corresponds 
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Figure 3: Solution of the dispersion relation of ( l96l) . The real part of N is plotted vs. k± 
for r = 1.11, u = 10.5, 6 = 0.10, s = —13.4, which corresponds to strong negative parallel 
velocity shear with 5 + s < 0. The energy signatures of the modes are indicated. 

to the drift-shear wave that has = iV+ at k± = 0. The second to upper is the Alfven 
wave, which has N = Nr at k± = corresponding to u/k = va, and the lowermost curve is 
its negative counterpart. The remaining wave is a drift-shear mode with = at k± = 0. 
For larger k± there exist regions of instability. 

For the stable degrees of freedom of Hamiltonian systems there exists a special form, a so- 
called normal form, to which all such systems can be mapped by a canonical transformation, 
{q,p) < — > {Q,P). The algorithm for this, which uses the real and imaginary parts of the 
linear eigenvectors, was first proven in total generality in [IH]. (See e.g. [22] for a more recent 
source for a version of the algorithm and [ST] where a plasma example is worked out.) For 
the stable modes this normal form Hamiltonian is given by 



where the prime on the sum means the ky and k^ values for the unstable modes are absent 
(recall the suppressed sum on k^) and the frequencies cu^*^ > 0. This Hamiltonian is merely 
that for a collection of simple harmonic oscillators, except for the presence of the signature 
a^^ G {1, —1}. Modes with signature a^^ = —1 oscillate, but are negative energy modes. 




(99) 



k 



i=l 
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The signatures of the modes can be obtained by inserting the eigenvectors into the Hamil- 
tonian (|88|) . If T^*'' is the eigenvector associated with the mode indexed by i and k, then 

its signature is given by the sign of Hl = T^*"* HlT^^ , where Hl is the matrix of the bihn- 
ear form (l88l) . To determine the signature of all modes, it is only necessary to do this for 
— oo and — 0. 

From (!97|) it follows for oo that uoK^/k ~ ±.va\J5{8 + + (5^). For these two 

modes in this limit, (!79|) is dominated by contributions from the terms | Jp, |Vv?p, f^, and 
D^, giving 

2 I 12 

i^L~^^^^, AT^O, K^^oo, (100) 

and thus they both have positive signature. Note, because s + 5 < 0, these modes are 
unstable for Fig. [3] and (11001) does not apply. For the remaining two modes the energy in 
this limit behaves as 

Hl- 2,2(m -^m^(n t/xr W AT ±iV, , ^ oo . (101) 

r^di{Ns =F N+){Ns =F A^-) 

From fllOip if follows that for the examples of Figs. [2] and [3l the mode approaching A^^ is 
negative while that approaching —Ns is positive. 

In the limit k,± the modes have the values ±Nr and N±, plus corrections to these 
values of order k\. For the Alfven waves that emerge from ±Nr the energy vanishes to leading 
order, reflecting the vanishing of line-bending in this long wavelength limit, but proceeding 
to order gives positive energies for all the stable Alfven waves of Figs. [2] and [31 Note, 
for Fig. [2], the signature of the negative Alfven wave — A^^ is consistent with the result of 
the K± ^ oo calculation: because this mode does not traverse zero or suffer a bifurcation 
to instability, it cannot change signature. Similarly we obtain that the drift-shear mode 
emerging from A^+ has negative energy for the values of both Figs. [2] and [21 while the A^_ 
mode has positive energy for Fig. [2l and negative energy for Fig. [31 As s goes from 18.5 to 
— 13.4, its frequency goes through zero and the mode changes signature. 

In Fig. [21 there is an unstable gap that lies between 2.39 < k,± < 4.91. At the value 
Kj_ ~ 2.39, we have a 'collision' where the frequencies of the Alfven wave and the upper 
drift-shear wave match and a transition to instability occurs. This transition is one of two 
types that occur in Hamiltonian systems; the other occurs at zero frequency as is the case 
for ideal MHD instabilities of static equilibria. A necessary condition for the transition to 
instability at nonzero frequency is that one of the modes must be a negative energy mode. 
This result is known as the Krein-Moser theorem [S2] . In Fig. [21 the transition that occurs 
at K± ~ 2.39 and its inverse that occurs at k± ~ 4.91 are both Krein-Moser transitions. 
Consequently one of the modes involved must be a negative energy mode, and we have 
shown this to be the case. It is important to note that the Hamiltonian formalism is the 
only reliable way to determine the existence of negative energy modes of a system. If one 
chooses an initial condition that only excites the drift-shear wave, the energy of the system 
will be negative. This is essentially what is shown when the sign of Hl = T^*"* Hi,Tf^ is 
obtained. 
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The existence of negative energy modes is a necessary condition for the transition to 
instabihty, but it is not a sufficient condition. Upon colhsion, modes may merely pass through 
each other and remain on the stable axis. The Hamiltonian is not a frame dependent quantity, 
and it is possible to change the signature of a mode with a frame change. For systems that 
possess a momentum invariant, the energy in the new frame is the same as the old, but 
with doppler shifted frequencies. Sometimes it is possible to remove a negative energy mode 
by this procedure and thus obtain a Liapunov type stability argument that precludes the 
transition to instability [531 El]- In any event, when a transition at finite frequency occurs, 
one is certain of the existence of a negative energy mode. 

After the transition, the energy drops to zero. Unstable modes always have zero en- 
ergy, and thus they have no signature. After the transition, the unstable gap modes have 
eigenvalues that occur in quartet form, ±u!r ± ujj, and the unstable normal form is 



where the double prime indicates a sum over /c- values of the gap, the indices 1 and 2 denote 
the upper Alfven and drift modes, and the indices 3 and 4 denote the remaining two stable 
modes. In the new coordinates, the total Hamiltonian is given by Hl = H'^ + H'^. Recall, in 
all of these sums, the sum over the x-components of the wave numbers have been suppressed. 

In Fig. O the unstable gap arising from the collision of the Alfven and drift modes occurs 
for values 2.12 < < 3.99, and then another collision occurs when a stable wave with 
positive signature collides with the lower drift-shear mode at k,± ~ 4.95. This lower mode 
changed from its positive energy value in Fig. [2] to a negative energy mode upon traversing 

= 0. For > 4.95, the normal form is like fll02p . except now the stable modes are the 
lower Alfven wave and the mode that limits to A'"^, and the other two modes are unstable 
for arbitrarily large k^. 

From Eqs. fllOOl) and (110 II) we see that if all modes are stable they will have positive 
energy if {Ng =]= N^){Ns =i= A^_) > for both signs. When this is the case, the system is 
energy stable, which can be traced back the the positive definiteness of 6^F. This kind 
of stability is sometimes called energy- Casimir stability (see [12] for references to the early 
plasma literature where this idea was introduced). Being energy stable means that our 
Hamiltonian can be brought into the form (1^^ for all k, with a^^^ > for all i and k. 

7 Collisionless conductivity and tearing modes 

In this section we consider another linear application. We show that the Jacobi identity 
of the bracket '[ , ]', an identity essential for the Hamiltonian description of Sec. [3] and its 
Jacobi identity for the bracket '{ , }', can be used to show that the azimuthal current density 
responds locally to the electric field on each fiux surface. That is, the current is proportional 




k 




(102) 



i=3,4 
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to the parallel electric field and is independent of its gradient across flux surfaces. This 
would not be true in the presence of particle diffusivity or in the presence of a finite electron 
gyroradius. The locality of the conduction is important as it allows the linearized system of 
equations to be reduced to a system that differs from that describing resistive tearing modes 
in a sheared slab by merely replacing the collisional conductivity a by a spatially dependent, 
AC coUisionless conductivity. 



7.1 CoUisionless conductivity 

We consider perturbations of a plasma slab in which all equilibrium fields vary only in the 
x-direction, i.e. the fields variables have the following form: 

ij = M^) + M^y'^'"'"'' + c.c. , (103) 

ip = -ux/k + ^k{x)e''''~'^y + C.C., (104) 

Z = azx + Zk{x)e''^'-'^y + c.c. , (105) 

V = a^x + Vk{x)e''''-'''y + c.c. , (106) 

where, to avoid clutter, below we drop the subscript k on the tilde variables. The above 
equations are the same as those of Sec. 16.21 except here we assume that ipo{x) of (11031) is 
a polynomial of at most quadratic order in x. This assumption makes it possible to apply 
the final result to two cases, namely the homogeneous equilibrium of Sec. 15.41 and the case 
of a magnetic equilibrium with a resonant surface at x = 0. The presence of an equilibrium 
flow in Eq. (11041) corresponds to choosing a reference frame moving with the phase velocity 
of the perturbation. This choice is convenient because it allows us to ignore the terms with 
the time derivatives in our calculations below. 
The Jacobi identity for [ , ] is 

b,[V',e]] + [e,b,^]] + [v^,[e,¥']] = o, (107) 

where ^ represents any of the field variable of fll03p - fll06p and consequently it has the form 
C = io{x) + ^{x) exp{iujt — iky) + c.c. Linearizing (I107P and retaining terms of first order 
gives 

where [f,g]L is a short-hand notation for the linearized Poisson bracket between two fields, 
e.g. [ip, = -ikiip'o^ - ^0^) exp{iujt - iky). 

By replacing with Z in fllOSp and using ([2]) one obtains 

o 

- I^IV-, Z]l + i^h + i^'^Q^W, V]l - cp[v, V]l) = 0, (109) 

where J = —'V'^ip. On the other hand by replacing ^ with v in fllOSp one gets 

^[.,^U = -^(a.|[^,^U + ^^|b,^]) =0. (110) 
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The latter expression can be used to replace d/dy[v,il)\L in fllU9l) and obtain 

(111) 

Making use of (jlj), (111 II) can be reformulated as 

-dp—\^,Z]L = dp ^ ^ . (112) 

Ohm's law ([1]), on the other hand, yields 

- d.^l^P, Z], = -^[^, i:], - dl^iy,, J]. (113) 

Given that in the linearized Poisson bracket the dependence on y enters only through the 
exponential, the derivative with respect to y amounts to a multiplication times —ik. Bearing 
this in mind and combining (11121) with (I113p . yields 



dpkuaz + dpipQa^cpk"^ — u;^ + ip'o^cik'^ 



(114) 



Now if one introduces in the linearized Poisson brackets the explicit expressions (I103p - (I106I) 
for the fields, one obtains 

ia{x) (cui! + ij'^k^) = -k'^i) + i)" (115) 

where 

uu^ + Cpdp^kuKH -uj'^ + i'o'^clk'^ , 

cr(x) := ^— o ^ — , 116 

^ ' iuj{k''%^dl + dl{-u^ + %\lk^)) ^ ' 

with uo^ = djskaz and ujkh = ka^. This shows that Ohm's law can be written as a propor- 
tionality relation between the amplitudes of the projection of the current density along z 
and of the electric field along the poloidal magnetic field. The quantity a{x) then plays the 
role of a spatially dependent conductivity. 

The case ipoix) = x'^/2Ls corresponds to an equilibrium with scale length Lg and with 
a resonant surface at a; = 0. If for this case we restrict to a thin layer around the resonant 
surface, the system comprised of (11151) and the vorticity equation ([3]) can be approximated 

by 

ia{x){uj^ + k\\ip) = i)" and uj^"{x) = , (117) 

where k\\ = kx/Lg and y-derivatives, being negligible in the layer, have been dropped. Thus 
the layer equations for the present model take the same form as those of MHD, 

unp"{x) = xip"{x) and cr(x)-E'|j = '?/'"(x) , (118) 
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except for the replacement of the conductivity by the spatially varying AC conductivity (t{x) 
of ffTTB]) . Here E\\ = i{uLitJj + x^p). 

The case ipoi^) = (^tp^ corresponds to the homogeneous equilibria of Sec. [61 Here a is 
constant and Eqs. flll7p become 

ia{ujtp + a^kip) = -k'^tjj + ijj" and uj{-k'^ip + ip") = -a^k{-k'^tjj + tjj"). (119) 

The solvability condition for this system, for solutions with dependence on x of the form 
exp(— ■jA;^;^), gives again the dispersion relation (!96|) . 



7.2 Collisionless tearing mode 

Now the dispersion relation derived above is used to obtain the growth rate for collisionless 
tearing modes. We restrict attention to the case of moderate A' where the constant-?/^ 
approximation applies. Mirnov et al.[55j have recently described the opposite case of large 
A' using a two- fluid model analogous to that of [36] . 

In the constant-?/' approximation, the dispersion relation for the tearing mode follows 
from the matching condition for the magnetic perturbation. 



oo 



A' = i / dxJ. (120) 



CO 



At the resonant surface, k\\ = 0, the conductivity is very high due to the high electron 
mobility, iuja{0) = {l—uj^/uj)/d1. Away from the resonant surface, however, the conductivity 
decreases rapidly due to the shielding of the electric fleld by the electron motion along 
the magnetic fleld. The shielding is described by the fcyc?^ term in the denominator. The 
region of high conductivity is called the current channel and for moderate tearing parameter 
A', it contains most of the current in the reconnection layer. In the current channel, the 
conductivity may be approximated by 

a(x)^-dj 



e 

iuj 



Substituting this in the matching integral and evaluating the integral gives 

A' = -in{uj - uj^) , ^\ , 
kdpde 

whence we obtain the dispersion relation 

UJ = UJ,+l — -^de. (121) 
ttLs 

In the limit /3 -C 1 (I121|) agrees with the kinetic result of [56], aside from a factor of 
2/y^ ~ 1.13 in the growth rate, and it agrees with the fluid result of ^57j in the low f3 limit. 
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8 Saturation of the collisionless tearing mode 



For our last application of the Hamiltonian formalism we use a Casimir invariant to find the 
nonlinear saturated state of the collisionless tearing mode following an approach similar to 
that in [58j, except that here we make use of the constant-'?/' approximation to simplify the 
analysis. For simplicity we consider the cold plasma limit where dp = cp = 0. Inspection 
of ([1]) reveals that the Casimir C2 of (125|) survives but with D replaced by tpe, Eq. ([3]) is 
unaltered, and Eqs. ([2]) and (HI) ensure that if initially Z, f = 0, then they will remain so. 

We use the invariance of C2 to describe what becomes of an unstable un-reconnected state, 
T^C') = Bqx^ /2Ls, such as that described in Sec. 17.21 as it evolves into a final approximate 
equilibrium state ip'^'^'^ = Bqx^ /2Ls + ^l) cosy that represents a magnetic island of half-width 
w = [ALgip / BqY^'^ . For convenience, throughout this section we use dimensionless units 
where lengths are scaled with w and the flux with Bq/Ls- In these units, for example, 
becomes 1})'^°^^ = x^/2 + cosy /A. 

Choosing ]C{ipe) = ^(V'e — '^e)/2vr singles out the surface of constant ipe, yielding 

where dxipeix^y^t) is to be evaluated at x = 4'e^ii'e,y,t), and this can be done at any 
time. From ip^^^ we obtain ip^^ = /2 — rf , where rf := dl/w"^] whence for the initial state 

^2(4) = 1/^2(4+^. 

Because the final state is an equilibrium state, ip^^"^ is a function of ip^°°^ according to 
ip^^ = ip^°°'^ + rf I [ip''°°^) ^ where / is the final current profile of the saturated island. Thus 
( 11221) becomes 

where dxip^°°\x,y) is to be evaluated at x = {ipi°°^)~^{ipe,y) and dipi°°^ / dip^°°\ being ex- 
pressible as a function of ipe alone, can be pulled outside the integral. The final current 
profile follows by setting the initial and final values of C2{ipe) equal at each ipe, yielding 

di^e 1 / dy ^^24) 



Although '?/'(°°) is unknown, for small A', dxip^^^ k, x ^ (4-?/'(°°) — cos y)/2. Integrating 
Eq. ([121) yields, 

4 = ,2(^(oo))_^2_ (^25) 

where 

i(7/>(~)) = — / dy a/4?/'(°°) - cosy. 
27r ./n 
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and where we have used the fact that hm^^oo(V'e('^) —'ip) = to set the integration constant 
to zero. 

Using = ip'^°°'> +'rfl{ip^°°^)^ and dropping the 'hat' and the label oo, gives an equation 
for the current profile, 

/(^) = + (126) 
The function is easily evaluated in terms of elliptic integrals: 



^ r 1,/W+WE{1/{2^ + 1/2)), for ^ > 1/4; 

\ (E(2^ + 1/2) + (2^ - l/2)ir(2^ + 1/2)) , for - 1/4 < ^ < 1/4 , 

where K and E are the complete elliptic integral of the first and second kind. 

Figure H] shows a comparison of the above current profile with that of Rutherford [38] for 
resistive diffusion. Observe, the profiles are qualitatively similar. Substituting the current 
profile in the matching relation of (11201) yields the saturation amplitude w = A'd^/G, where 

, , ^, , , f dy cos-u 
dipli^j) (p = 0.19 



G = 8 

This is close to the value G = 0.205 obtained by Drake and Lee [59] using a kinetic model. 
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Figure 4: Comparison of the current density profiles, one determined by parallel electron 
momentum conservation and the other by Ohmic diffusion (dashed line) 



9 Summary and conclusions 

In the early sections of this paper we presented the noncanonical Hamiltonian formulation 
of the four- field model of [36], and showed that the associated Lie-Poisson bracket has four 
new independent families of Casimir invariants. These invariants led us to the discovery of 
variables in which the Poisson bracket has the simple form and in which the system 
can be written in the compact form of (1351) - (1371) . 
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In Sec. O we used the Hamiltonian formulation to obtain a variational principle that 
gives a set of coupled differential equations that generalize the Grad-Shafranov equilibrium 
equation. In the limit of vanishing electron mass {de —>■ 0) the equilibrium equations reduce 
to previously known results. This limit provides some insight into the relationship of the 
Casimirs to the more familiar conserved quantities of conventional low-/? drift models. We 
have presented two solutions of the equilibrium equations, the first describing dipole-like 
equilibria and the second describing homogeneous equilibria that support drift-acoustic and 
Alfven modes. 

In Sec. |6] we investigated the linear dispersion relation for homogeneous equilibria and 
described the map to the appropriate Hamiltonian forms. We also presented thresholds for 
spectral and energy stability. In Sec. 16.31 we described a method for determining the energy 
signature of a mode. This method is of general utility and can be applied to all valid models, 
provided one understands their Hamiltonian structure. In obtaining reduced fluid models, 
there can be ambiguity about the energy for the full dynamics, and linear theory alone cannot 
be used to uniquely determine the correct energy of the linear dynamics. The only reliable 
way to determine the energy is from a Hamiltonian or action principle formulation, where 
the energy for the linear dynamics is obtained by expansion of a Hamiltonian associated 
with time translation symmetry. 

In Sees. [7] and [8] we demonstrated the usefulness of the Hamiltonian formulation for the 
analysis of the linear collisionless tearing mode and its nonlinear saturation. In the case of 
linear stability, the Jacobi identity allows the reduction of the system to a form analogous to 
that of MHD but where the conductivity is replaced by a spatially varying AC conductivity. 
Applications of the formalism left for future work include the study of the stability of the 
saturated states against secondary modes and an investigation of saturation in a more general 
dynamical context using additional Casimirs. 

Another area for future work concerns the families associated with the invariants T±, 
which generalize a pair of Casimirs G± found for a low-/? two-field model derived in [1]. A 
natural question is whether the invariants T± play a role analogous to the one played by 
G± in the two-field limit in determining the alignment of current density and vorticity along 
the separatrices of the magnetic field during the nonlinear evolution of the system [31] . The 
present model makes possible an investigation of this question along the lines carried out in 
[22]. 
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